function p = newton_raphson_reciprocal(D,N)
    p = 1;
    Dp = D;
    for i=1:N,
        p = p*(2-Dp);
        Dp = Dp*(2-Dp);
    end
end

